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Abstract. High resolution, mid-infrared (MIR) images of a set of nine nearby active galaxies are presented. The data were 
obtained with the TIMMI2 instrument mounted at the ESO 3.6 m telescope using a set of N-band narrow filters. The resulting 
images have an angular resolution of 0.6"-l". The MIR emission has been resolved in four galaxies: NGC253, NGC1365, 
NGC1808 and NGC7469. The images unveil a circumnuclear population of unknown MIR sources in NGC1365 and NGC1808, 
coincident with radio sources. These MIR/radio sources are interpreted in terms of embedded young star clusters. A high- 
resolution MIR map of NGC253 is also presented, and enables the identification of a previously unknown MIR counterpart to 
the radio nucleus. Extended MIR emission is detected in NGC7469, and concurs with previous observations in the NIR and 
radio. For this source, an interesting morphological difference between the 10.4 jim and the 11.9 yum emission is observed, 
suggesting the presence of a dust-rich micro-bar. Our MIR images of Circinus do not show resolved emission from the nucleus 
down to an angular scale of 0.5". In the case of NGC2992, an upper limit to the extended MIR emission can be set. Finally, we 
provide new MIR flux measurements for the unresolved AGN in NGC5995, IZwl and IIZwl36. 

Key words. Galaxies: individual: NGC253, NGC1808, NGC1365; Galaxies: star clusters; Infrared: galaxies; Galaxies: Seyfert 
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1. Introduction 

As mid-infrared (MIR) radiation is much less affected by ex- 
tinction than visible light, dust embedded objects in the nuclear 
regions of galaxies are more efficiently detected in this wave- 
length range. High spatial resolution in the MIR requires the 
use of large-aperture telescopes and thus is restricted to ground- 
based instruments that provide a substantial gain in spatial res- 
olution over currently available MIR instruments in space. The 
gain in spatial resolution is, in the context of active galaxies, 
essential since in many cases the relative contributions to their 
total power of the starburst activity and of the active galactic 
nucleus (AGN) are unknown. Let us briefly discuss how high 
resolution MIR observations can help probe/disentangle these 
two phenomena. 

Young bursts of star formation are likely to be substan- 
tially obscured by dust. Therefore, high resolution MIR imag- 
ing is an appropriate means of identifying them. Indeed, re- 
cent infrared (IR) observations have revealed the existence of 
very bright embedded star clusters, containing few thousands 
O stars, which are invisible on images obtained in the optical. 
The first object of this type to be discovered wa s the deep! 
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buried star cluster in the Antennae galaxy by Mirabel et 

It produces 20% of the total MIR emission of the 
galaxy, and contains 1.6 x 1O 7 M of stars hidden behind an 



Send offprint requests to: E. Galliano, e-mail: egallian@eso.org 



A v = 10 cloud of dust JGilbert et alJl200fj). S imilar objects 
have been found in NGC5253 dGoriian et alfeOOll SBS 033 5- 
052 frlante & SaiivagefeOOA and II7w40 feeck et alfcOOA. 

High resolution MIR imaging is also essential for the iden- 
tification and classification of AGN. For sources affected by in- 
trinsically large extinction, an optic ally-based classificati on can 
be quite misleading. For example, Krabb e et alJ (1200 II) noted 
that the classification of "active galaxies" based solely upon 
their optical properties is a poor indicator since X-ray or radio 
observations could reveal the presence of an AGN in "classi- 
cal" s tarburst galaxies. Conversely, ISO observations jLaurenll 
Il999h of optically classified AGN show that their MIR emis- 
sion is often dominated by star forming regions rather than 
by the AGN itself. Thus by spatially resolving emission from 
star forming regions and AGN, and photometrically measuring 
each through high resolution imaging, one can determine their 
relative contributions to the total integrated power. For AGN, 
precise and spatially resolved MIR photometry is mandatory 
in order to construct spectral energy distributions (SED) that 
are free from starburst contributions. The resultant SED may 
then be compared to models developed in the framework of 
the unified scheme. Due to the intrinsic small size of the dusty 
torus in AGN (1 to a few xlOOpc), only a few can be di- 
rectly probed through a spatially resolved mapping of their 
MIR emission. One of these is the extensively studied type 2 
AGN in NGC1068. MIR imaging of this source at 0.6" reso- 
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lution revealed extended and collimated 10 and 20 //m emis- 
sion up to 5" from th e core, along the direction of the radio 
jet jAlloin et al.ll2.Qpoj) . At higher angular resolution (0.1 " af- 
ter deconvolution: iBock et aljfeOOOtlTomono et al.ll200lh . the 
complex structure of the MIR emission is foun d to be coinci- 
dent w ith the narrow line region (NLR) clouds dGalliano et alJ 
120031) . Therefore, the access to high resolution imaging brings 
new clues regarding the morphology of the material around the 
central engines of AGN. 

This paper presents high angular resolution MIR images 
of a sample of active galaxies, obtained with the TIMMI 2 in- 
strument at the ESO 3.6 m telescope, La Silla observatory. 
Observations and data reduction are presented and discussed 
in Sect. [2] and the results are given in Sect. [3] First, the gen- 
eral results such as flux measurements and maps for the re- 
solved objects are given in Sect. 13.11 Then, we provide more 
detailed analyses for NGC253 (Sect.ED. Circinus (Sect.E), 
NGC1808 (Sect. E3), NGC1365 (Sect. and NGC7469 
(Sect. l3~6t . Sect.|4]gathers high spatial resolution IR flux mea- 
surements and provides SED for the targets. Finally, we discuss 
in Sect. [5]the newly discovered populations of MIR sources in 
NGC1808 and NGC1365: a simple model supports the inter- 
pretation in terms of embedded star clusters. A detailed discus- 
sion of these sources will be presented along with new obser- 
vations in a forthcoming paper. Summary and conclusions are 
given in Sect[6] A Hubble constant Ho = 75 kms -1 Mpc~' is 
used throughout the article to compute the distance to the ob- 
jects, except when measured with a more precise method. 

2. Observations, data reduction and templates 

2.1. Observations 

The original aims of this study were two-fold: (1) perform pre- 
liminary observations of active galaxies in order to prepare 
future observations with VISIR at the VLT, (2) look for ex- 
tended MIR emission directly related to the AGN , as was ob- 
served along the radio jet direction in NGC1068 (lAlloin et alJ 
2000: IBock et al l2000l iTomono et all200lh . To address these 
aims, we chose a sample of nearby (few Mpc) to moderately 
nearby (few hundreds of Mpc) AGN. The observations quickly 
showed interesting by-products such as the star forming envi- 
ronment around some AGN: the MIR counterparts of known 
radio sources in NGC1365 and NGC1808 became unveiled. 
After these first findings, we decided to include in our sam- 
ple NGC253, the archetypal starburst galaxy. Table [2 displays 
basic information about the targets. 

The observations were performed with TIMMI 2 (Thermal 
Infrared Multi Mode Instrument 1 ) installed at the 3.6m tele- 
scope on La Silla during four runs, on 2001 March 10-11, July 
5, October 8 and November 26. TIMMI 2 is equipped with a 
240 x 320 pixel AsSi BIB detector and we observed with an 
0.2" pixel size, resulting in a field of view of 48"x 64". 

Four narrow filters spanning the N-band were avail- 
able: the 8.9 /an filter (7.90-9.46 /mi), the 10.4 /mi filter 
(9.46-1 1.21 jum), the 11.9 //m filter (10.61-12.50 //m) and the 

1 http://www.ls.eso.org/lasilla/Telescopes/360cat/timmi/ 



12.9 /mi filter (1 1.54-12.98 /mi). The labels of the filters do not 
correspond to their effective wavelengths, but rather refer to the 
spectral features they probe. One filter is also available in the 
Q-band (20 //m). The seeing was good, ~0.6-1.0" in the visible 
and allowed to obtain diffraction-limited images at all observed 
MIR wavelengths. The main difficulties in performing MIR ob- 
servations from the ground are the high level of sky emission 
and the variability of the sky transparency, particularly in the 
20 //m window. 

2.2. Data reduction 

The standard chopping/nodding technique was used to deter- 
mine the sky and telescope background for subtraction. We ob- 
served in "small imaging mode", which means that the nodding 
direction was perpendicular to the chopping direction. On-line 
reduction at the telescope allowed the chopping throws and 
nodding offsets to be adjusted depending on the extension of 
the sources. Chopping throws between 15" and 20" and nod- 
ding offsets from 15" to 30" were used. MIR standard stars 
were observed before and after each source observation in 
order to calibrate the photometry of the final images, evalu- 
ate variations of the sky transparency, and estimate the Point 
Spread Functions (PSF) for further image deconvolution. IDL 
routines were written and used for data reduction. Glitches 
were carefully removed and the data were corrected for the 
TIMMI 2 bias: each of the sixteen readout channels of the chip 
has a response which varies with time. This appears on the im- 
ages as sixteen vertical bands with slightly different intensity 
levels. This detector feature makes the flat field correction use- 
less, as it assumes a stable detector response, at least over the 
time scale of an observing night. These reduction steps were 
performed on each plane of the data-cubes, each of these planes 
corresponding to the subtraction of two consecutively chopped 
images. Eventually, the planes of each reduced data-cube were 
co-added resulting in the scientifically exploitable image. 

Where resolved structures were present, we have de- 
convolved the images in order to better separate the dif- 
ferent sources. Deconvolution was performed using the 
mr deconv routine of the M R/1 software package developed 
by Murtagh & Starck ( 1999). Standard stars were used to deter- 
mine the PSF. The MR/ 1 software package performs an efficient 
regularisation in the wavelet space. This procedure is excellent 
at pulling low signal out of noise, clarifying the distribution of 
sources. The results obtained after deconvolution were suitable 
for comparison with high resolution cm radio images and we 
could indeed clearly identify MIR counterparts for a number of 
radio sources. For low brightness sources within the resolved 
structures, we measured the fluxes on the deconvolved images. 
The precision of this procedure was tested with the following 
simulation: an artificial MIR image of NGC1808 was created, 
in which weak point sources were given the fluxes measured 
on the deconvolved image. Poisson noise corresponding to the 
high background level was added. Then, this fake image was 
deconvolved: after deconvolution, the input fluxes were recov- 
ered with an accuracy better than 10%. The brightest among the 
weak sources were all recovered, while a few false very weak 
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Table 1. General properties of the sample galaxies: coordinates, distance, AGN type, IR flux 



Galaxy 


RA 


Decl. 


Dist. 


Scale 


Activity 


IRAS 12pmflux c 




(J2000) 


(J2000) 


Mpc 


pc/arcsec 




Jy 


NGC253" 


00 h 47 m 33.1 s 


-25° 17' 18" 


2.5 


12 


Starburst + AGN? 


24.0 


Circinus 


14 h 13 m 09.3 s 


-65°20'21" 


4.0 


19 


Type 2 AGN 


18.8 


NGC1808* 


05 h 07 m 42.3 s 


-37°30'46" 


10.9 


53 


Starburst + AGN? 


4.4 


NGC1365 


03 h 33 m 36.4 s 


-36°08'25" 


18.6 


90 


Type 2 AGN 


3.4 


NGC2992 


21 h 32 m 27.8 s 


+ 10°08'19" 


30.1 


150 


Type 2 AGN 




NGC7469 


23 h 03 m 15.6 s 


+08°52'26" 


66 


320 


Type 1 AGN 


1.35 


NGC5995 


15 h 48 m 24.9 s 


-13°45'28" 


100 


490 


Type 2 AGN 


0.39 


IZwl 


00 h 53 m 34.9 s 


+ 12°41'36" 


245 


1186 


Type 1 AGN 


0.51 


IIZwl36 


21 h 32 m 27.8 s 


+ 10°08'19" 


252 


1222 


Type 1 AGN 


<0.2 



" A discussion about the nature of the radio nucleus of NGC253 can be found in lForbesetallj200d). 

* The AGN nature of the nucleus of NGC1 808 is not cl early established, see the discussion in Krabbe et alJfeOOll) and i 
c Flux of the entire galaxy, from lMoshir & et all fi 990). 
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sources appeared on the deconvolved image, but could be eas- 
ily identified and excluded since they were more compact than 
real sources. 

2.3. Photometric error estimation 

The uncertainty on the flux measurements is dominated by vari- 
ations of the sky transparency. In the 10 /mi window, variations 
of the photometric calibration (ADU to flux) during the night 
were found to be of about 10%. One extreme case of trans- 
parency variation occured on 2001 November 26: the photo- 
metric calibration varied by 50% in about 20 min probably be- 
cause of spatial variations of the sky transparency. This hap- 
pened during the observing sequence for NGC1365 through 
the 12.9//m filter. The corresponding data-cube was checked: 
we split it into equal slices, for each of which the source fluxes 
were measured and found to be equal. This ensured that the 
sky transparency had not changed dramatically while perform- 
ing the science observation, but rather differed between the two 
standard star observations. Unfortunately, there is no possibil- 
ity to check the spatial variation of the sky transparency. In 
conclusion, we will keep in mind that the 12.9 /mi flux mea- 
sured for NGC1365 is uncertain. One must also consider that 
the values for the standard star fluxes are known to a precision 
that is not better than 10%. Finally, we adopt a formal error of 
+20% on the quoted flux measurements. 

2.4. Templates for MIR sources: AGN, HII region and 
PDR 

The MIR emission observed in galaxies has different origins: 

- The Rayleigh-Jeans emission of the evolved stellar popu- 
lation is the dominant contribution in early-type galaxies 
which are poor in gas and dust. 

- The dust, irradiated and heated by an underlying contin- 
uum, emits a MIR continuum. Large grains are in thermal 
equilibrium with the radiation field and have a blackbody- 
like emission. There also exists a population of very small 
grains (size < 50A) which never reach the thermal equilib- 
rium. These very small grains (hereafter VSG) can get to 



high temperatures (T > 100K) after absorbing one single 
UV photon. The grains cool down by emission longwards 
of 10//m before they absorb new photons. The VSG emis- 
sion is typical of regions which are actively forming stars, 
such as HII regions. 

- The ionised gas of the interstellar medium emits several 
strong forbidden lines in the MIR. Lines from Ne, Ar, Si 
as well as S and Fe are commonly observed in the spec- 
tra of AGN (excluding BLLac and OVV) and HII regions. 
In particular, [Arlll] 8.9/rni, [SilV] 10.4 /mi and [Nell] 
12.8 /mi lie within the wavelength range of the set of fil- 
ters we used for this study. 

- The unidentified infrared bands (UIB) observed at 6.2, 
7.7, 8.6, 11.3 and 12.7/mi are characteristic of photo- 
dissociation regions (hereafter PDR), where molecules 
form and their physical conditions depe nd on the UV 
photo n density. Since the original work of iLeger & Pugetl 
( 1984§, these bands are interpreted as the bending and 
stretching vibration modes of the C=C and C-H links in 
polycyclic aromatic hydrocarbons (PAH). 

In Fig. [H template spectra over the 7-14/im range 
dLaurent et alJl200ol) are displayed for the three main types of 
MIR sources: AGN, HII region and PDR. These templates are 
useful for comparison with observations: we therefore over- 
plotted on the template spectra the flux densities that would 
be measured through the TIMMI 2 filters. 

The AGN spectrum is dominated by continuum emis- 
sion from the dusty torus, and hard silicate absorption around 
9.7 //m can be observed. Usually, PAH emission is absent or 
very weak since the carriers are destroyed by the strong radia- 
tion field. 

The HII region spectrum displays a steeply rising con- 
tinuum, due to VSG emission. The ionised gas produces the 
strong [Arlll] 8.9 /mi, [SilV] 10.4/miand [Nell] 12.8 /mi, but 
the PAH carriers are destroyed by the radiation field as in the 
case of AGN. 

The PDR spectrum is dominated by the emission of the 
PAH. The PDR forms at a distance from the hot stars which is 
large enough for the UV continuum to be weakened and soft- 
ened; this allows molecules to be formed in the cool gas. 
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AGN template 



"v. 



7 8 9 10 11 12 13 14 

wavelength [fim] 

HII region template 




7 8 9 10 11 12 13 14 

wavelength [/xm] 

PDR template 




7 8 9 10 11 12 13 14 

wavelength [/tm] 



Fig. 1. Template spectra for the 7- 14 /mi emission of the three 
main types of MI R sources: AGN, H II region and PDR. The 
spectra are from Lau rent et alJ 11200(1) . The corresponding flux 
densities through the 8. 9 //m, 10.4 /mi, 11. 9 /mi and 12.9 /mi fil- 
ters are overplotted: the horizontal bar marks the width of the 
filter, while the abscissa of the dots correspond to the wave- 
lengths of the TIMMI2 filter names. 



We shall use the flux ratios computed with these templates 
to distinguish between the various object types and to help in- 
terpret the MIR emission from our sample sources. 

3. Results 

An overview of the observational results is presented in this 
section. In the first part, we give results for the whole sam- 
ple, and provide a short general discussion. Then we discuss 
in more detail the most interesting cases. Table |2] summarises 



the full set of measurements. For each target, we provide the 
fluxes and relative positions for all the identified sources. The 
targets are listed by increasing distance. MIR extended emis- 
sion (resolved and/or diffuse) is detected in 4 objects: NGC253, 
NGC1808, NGC1365 and NGC7469. Contour maps of the 
emission for these sources are presented in Fig. [2 For each ob- 
ject, the brightest source (which is the AGN except in the case 
of NGC253) is labelled Ml, while the circumnuclear sources 
are labelled M2, M3, etc... 

3.1. Overview 

MIR emission was detected for all the targets in the sample: 

- In NGC253, six distinct sources have been identified in the 
central region. This provides a gain in spatial resolution 
with respect to previous observations, which helps clarify 
the distribution of the MIR emission from this galaxy. 

- In Circinus, emission has been detected in all four narrow 
N-band filters. This is also the only target that we have de- 
tected at 20 //m. 

- In NGC1808 and NGC1365, a previously unknown popu- 
lation of MIR sources in the inner 10" region around their 
AGN is unveiled. 

- In NGC7469, extended diffuse emission is detected around 
the AGN and shows interesting morphological peculiari- 
ties. While the 12.9 //m map shows an emission quite sym- 
metrically distributed around the nucleus, the 10.4 /mi map 
unveils the presence of a bar-like structure along the N-S 
direction. 

- The last four objects, NGC1992, NGC5995, IZwl and 
IIZwl36 are detected, each through two to three filters, as 
single unresolved sources. 

The measured fluxes are generally in good agreement 
with previous measurements (see dedicated sections hereafter). 
Most previous measurements had been derived from data ob- 
tained with poorer spatial resolution and from photometry per- 
formed through much larger apertures than ours, hence encom- 
passing a much larger fraction of the galaxies than the nuclear 
regions we are probing with TIMMI2. Therefore, this suggests 
that most of the MIR emission from these targets is indeed aris- 
ing in their central ~ 1 0" region and can be accounted for by the 
AGN and the discrete sources revealed in the TIMMI2images. 

In none of the targets do we find extended MIR emission 
compara ble to that observed in NGC1068 or 3C 66B along the 
radio jet ( Allo in et alJbOOOtlTanslev et all2 000). except a ten- 
tative detection in NGC1365 (see Sect. 13. 5i . According to the 
AGN unified scheme, we were expecting to detect collimated 
MIR emission in the type 2 objects of our sample. It is not so. 
Still, we believe that an analysis on a larger sample and from 
images with higher sensitivity is mandatory to understand the 
diversity of MIR morphologies of type 2 objects, and is left for 
future study. 

In the following, we present, in some detail, the re- 
sults obtained for NGC253, Circinus, NGC1808, NGC1365, 
NGC7469. 
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NGC353, 11.9pm 



NGC7469, 10.4pm 



NGC7469, 11.9/im 
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Fig. 2. Deconvolved maps for the objects with resolved emission. North is up and East left. Dotted contours correspond to lcr of 
the background. Source labelling is according to the scheme described in Sect. [5] Galaxy and filter identifications are given as 
a title for each panel. The lower left square on each panel displays the corresponding PSF shape (FWHM in the E-W and N-S 
directions). 



3.2. NGC253: MIR detection of the radio nucleus? 

NGC253 hosts a vigorous starburst in its inne r lOOpc cen- 
tral re gion and is the archetypal starburst galaxy. iTurner & Hoi 
( 1985) discovered a set of highl y aligned compact radio source s 
along the -NE-SW direction, luivestad & Antonuccil Jl997h . 
performed high resolution and sensitive VLA observations 
and identified 64 sources. These sources probably belong to 



a nuclear ring seen at high inclination, and thus appear in 
a string-like distribution. The most powerful of these radio 
sources (TH-2) is considered to b e the true nucleus of NGC253, 
and is thought to be an AGN llForbes et al. 2000). Imaging 



of NGC253 


has been 


carried out in the NIR ( Forbes et alJ 


1991, 


1993; Sams etal. 


Il994t iKalas & Wvnn-Williamslll994 


and MIR (Pina et alll992llKeto et alll993llB6ker et alJll998: 



These observations show that the NIR/MIR 
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Fig. 3. NGC253, two tentative registrations of the MIR 
map (contours) with respect to the radio map (brightest 
sour ces shown by circle s ). The radio map was obtained 
by llJlvestad & AntonuccH dl997l) at 2 cm. (a) Absolute as- 
trometry of the NIR/MIR images. The three boxes show 
the absolute (B1950) positions and positional errors for the 
bright NIR/MIR source Ml: the dashed, solid and dotted 
line boxes show the respective regi strat ions of [P ina et al. 
dl992l) . lKalas & Wvnn-Willilmisld 19941) and lKeto et alJdl993l) . 
The small squares (0.2" side) record the positions of the 
TIMMI2 sources with the brightest source Ml assigned at 
the most probab le absolute position of peakl as derived by 
Keto et al. I dl993l) . The cross marks the absolute position of the 
radio nucleus, (b) Blind registration: Ml is positioned at the 
coordinates derived bv lKeto et alJ dl993l) . (c) Best registration: 
tentative alignment of the MIR and radio peaks. As in Fig. [2] 
the lower left square on each panel displays the PSF. 



emission consists of a bright, barely elongated source plus a 
secondary, fainter source to the NE of the brightest one. The 
brightest NIR/MIR source is believed to be a young massive 
star cluster llKeto et alll 999). Both sources are inside an elon- 
gated envelope, with PA similar to that of the radio string. 
However, a general lack of correlation between the IR and the 
radio sources is observed. 

The TIMMI2 image of NGC253 (through the 1 1.9 jum fil- 
ter) is the highest resolution MIR image of this galaxy available 
to date. The MIR morphology revealed in this image is fully 
consistent with the previous MIR observations at lower angular 
resolution, but the distribution of the emission is resolved into 
six discrete sources. The flux measurement at 1 1 .7 /mi made 
by IPina et alJ ll 19921) of Ml is best suited for comparison 
with our 11.9jum flux. Their flux of 4.8Jy is in good agree- 
ment with ours (5.1+0. 1 Jy). Comparison cannot be made with 
fluxes in the other filters bec ause of the presence of strong 
emission lines and PAH bands dKalas & Wvnn- Williams! 1 994 
iKetoet all! 9991). 

Again, at the resolution achieved with TIMMI 2, the source 
distributions in the radio and in the MIR are not obviously co- 
incident, even though the global position angle of both distri- 
butions is the same. 

We propose two registra tions of our M I R m ap with respect 
to the 2 cm map bv llJlvestad & Antonuccil dl997): 

(a) We use th e absolute registration of Ml performed by 
iKeto et al. (1993). It is consistent with other measurements 
from lPina etalJ i 19921) and iKalas & Wvnn- Williams! dl994) . 
and has the smallest uncertainty (about 0.4"): the measurement 
of the position of Ml gives a(1950)=00 h 45 m 5.630 s + 0.026 s & 
£(1950)=-25°33'40.57"+0.44". The resulting superposition of 
the MIR map and the radio sources is displayed in Fig. |3j?. We 
do not believe this blind registration to be the correct one, since 
the MIR string and the radio string are shifted one with respect 
to the other. However, we note that under this registration, M4 
is coincident with the radio nucleus. 

(b) We try to align the MIR and radio peaks. To do so, 
we simply minimise the sum, over all the MIR sources, of 
the distances to the closest radio source. Fig.|5J; shows the re- 
sulting map. This registration suggests, this time, that source 
M3 is the counterpart of the radio nucleus. The superposi- 
tion is globally much more convincing than the blind one. 
In this case, the B1950 coordinates of Ml can be derived: 
a=00 h 45 m 5.703 s & <5=-25°33'40.57". These coordinates are 
not consistent, within the qu oted error boxes, w ith the absolute 
coordinates of Ml derived by Ket o et alJ 111 9931). but are co nsis- 
tent with those, with larger erro r bar, bv lPina et alJ (fl992) and 
IKalas & Wvnn-Williamsldl994l) . 

Despite the inconsistency with Keto's measurement, we 
favour the second registration, for which the string of radio 
sources is well aligned with the MIR emission. In none of these 
registrations does source Ml have a radio counter part. The lack 
of rad io emission from Ml, already discussed bv lForbes et alJ 
(2000), is confirmed by our study and remains mysterious. And 
finally, our preferred registration places the radio nucleus in co- 
incidence with the MIR source M3. 
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3.3. Circinus: still point like at 0.6" resolution 



Circinus and HR5459 8.9/j.m profiles 
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Fig. 4. Photometric measurements for the Circinus galaxy. The 
fille d circles show our measurements. The ISO spectrum from 
LSturm et"al] J200 0) is overplotted for comparison. 



Circinus is the closest type 2 AGN, at a distance of only 
4Mpc. The nuclear region harbours a hea vily obscured A GN 
(N H > 10 24 cnr 2 based on X-ray data bv lMatt et al1ll999h . It 
also exhibits a prominent ionisation cone extending to the NW, 
and a starbur st ring at a radius of about 10" (~100pc) from 
the nucleus (tDavies et al . 1998). Circinus has been previously 
imag ed at high resoluti on with TIMM I ( Siebe nmorgen et all 
1 1997b and MANIAC dKrabbe et alJ kOOll) . both showing 
marginal extension of the central source. The latter study also 
revealed the presence of a low surface brightness diffuse halo 
of radius ~ 10". 

Aperture photometry was performed to obtain the N-band 
fluxes for Circinus within an aperture of 2" radius. This ra- 
dius was optimised using a curve of growth analysis. The 
background was measured in an annulus around this aper- 
ture (with inner and outer radii 2" and 4", respectively). The 
TIMMI 2 data points are plo tted in Fig. HI The 14" x20" aper- 
ture ISO-SWS spectrum bv ISturmet all (12000) is shown.The 
N-band TIMMI 2 measurements are in good agreement with 
the ISO spectrum, showing that the entire MIR emission comes 
from the nuclear source. At 20 yum, the discrepancy between the 
TIMMI 2 data and the ISO spectrum is difficult to interpret in a 
decisive manner. Observing and calibrating ground based data 
in this window is a difficult task. Moreover, the 20 yum emis- 
sion is expected to be more extended than the 10/i/m one, since 
it traces cooler material. Undetected diffuse emission on a scale 
of the order of 10", could be responsible for this difference. 

On the TIMMI 2 data, the comparison between the AGN 
and the corresponding standard star profiles does not reveal 
any sign of resolved emission in Circinus on a scale 0.3- 
0.4" (5.5-7.5 pc). The TIMMI2 image reaches a Icr back- 
ground level of 25 mJy arcsec~ 2 , which is insufficient to detect 
the broad N-band halo emission claimed by Krabbe et al. in 
2001 (17mJyarcsec- 2 ). 
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Circinus and HR5459 12.9/^m profiles 




Fig. 5. Circinus: profiles along the E-W direction at 
8.9 fim (top) and 12.9 jjm (bottom). The squares feature the pro- 
file of Circinus and its sampling, while the dashed and dotted 
lines show the profiles of the standard star HR5459, acquired 
before and after Circinus, respectively. 



3.4. NGC1808: MIR detection of radio sources 

NGC1808 is a morphologically pecu liar spiral galaxy (Sb c 
pec) located at a distance of 10.9 Mpc (Koribal ski et all 19 96). 
The n ature of its nuclear source i s not yet clearly estab- 
lished. IVeron-Cettv & Yfo onl jl9^h concl uded to the pres- 
ence of an AGN, while Forbes et al. il 19921) proposed instead 
that the nucleus consists of a SNR embedded within a metal- 
rich HII region. NGC1808 has been detected in hard X- 
rays ( Awa ki & Kovamall 19931 lJunkes et alJI 19951) . suggesting 
that an embedd ed AGN is indeed p resent in the nuclear re- 
gion. However, iKrabbe et al.l J200ll) show that, in a MIR vs. 
X-ray flux diagram, the location of NGC1808 argues against 
the presence of an AGN. 

The central region (inner ~750pc) of NGC1808 is under- 
going an episode of enhanced star formation. There were early 
reports on the presence of "hot spots" in its complex nuclear re- 
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NGC1808: radio, MIR and optical 
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Fig. 6. NGC1808: from left to right, radio, MIR and optical images. The map at 3.6 cm is from ICollison et all Jl994l) : our 
TIMMI2 map at 12.9 yum is displayed with a non-linear flux scale in order to highlight faint structures; the HST image (F658N) 
is from the HST archive. On the MIR and optical maps, the locations of the main radio sources are marked by circles of radius 
0.8". 



gion by Morgan ( 1958) and lSersic & Pastorizal d 1965b. Strong 
PAH emission was observed bv ISiebenmorgen et alJ J2001 ) 
from the star-forming region. Radio c ontinuum observations 
( Saikia et alJl99(A |Collison et al. 1994) revealed a population 
of compact radio sources, most of w hich do not coincid e with 
the optical hot spots. According to Saiki a et all |l990), they 
are more likely related to young SNRs than due to the thermal 
emission from giant H IT regions. A high resolution K -band im- 
age was published by iTacconi-Garman et alJ (1996), showing 
the presence of many distinct and compact (< 0.7") sources 
spread all over the circumnuclear region. These authors ar- 
gue that such sources are likely to be clusters containing hun- 
dreds of massive red super-giants. Not all of the radio emi tting 
source s have counterparts in the map bv lTacconi-Garman et alJ 
Jl996l) . Particularly, the three brightest circumnuclear radio 
sources have no NIR counterpart. 

Krab beetai1(l200lh found a nuclear flux of 0.33 Jy through 
an N-band (8-13//m) filter. This value is in reasonable agree- 
ment with ours (see Table although a precise compari- 
son is not possible because of the different filters used. The 
12.9/im map of the central region of NGC1808 is repro- 
duced in Fig. |6l together with the 3.6cm image obtained by 
Collis on et al . (1994) and an HST archive image in the opti- 
cal (F658N). The radio and MIR maps are strikingly similar. 
The most intense radio sources are featured by circles of ar- 
bitrary 0.8" radius on the MIR and optical maps. Several of 
the intense radio sources have a counterpart on the MIR image. 
Particularly, this comparison shows that M2, M3, M6 and M8 
are also radio sources. On the HST image, counterparts for the 
MIR/radio sources are not clearly detected, but may be identi- 
fied to faint sources that can be seen inside the circles in Fig. [6] 
M8 is the source that has most likely an optical counterpart. A 
comparison between a higher resolution MIR map and the HST 
image is needed in order to investigate which optical sources 
are the counterparts of the MIR/radio sources. Still, it is clear 
that the regions of most intense optical emission (the hot spots) 



do not coincide with either radio or MIR sources. The nature of 
these MIR/radio sources will be discussed in detail in Sect. [5] 



3.5. NGC1365: MIR detection of radio sources 

A global picture of NGC1 365 can be found in the dedicated 
review by Lindblad (1999). This spiral galaxy with a promi- 
nent b ar has been classified as SBb(s)I dSandage & Tammann 
Il98ll). An estimate of its distance is discussed in iLind blad 
( 1999J) and is currently set to 18.6+0.6 Mpc. At this distance, 
1" corresponds to 90 pc. The nuclear region harbours an AGN 
witnessed by the presence of both broad and narrow Balmer 
emission lines, as well a s narrow and strong forbidden e missio n 
lines of high excitation dVeron et alJll980tlAlloin et alJll98ll) . 
It is classified as a type 1 .5 AGN, but varies from type 1 to type 
2 according to the level of its hard radiation emission. 

It has long bee n known to show interest ing structures in 
its nuclear region dSersic & Pastorizal 1 19651) . In the optical, 
this region shows several bright hot spots within 1 kpc around 
the AGN. The HST images resolve the circumnuclear re- 
gion into a large numbe r of bright, compact star clusters o f 
only a few parsec radii jRristen et alJll997t lLindbl"adl [l"999 ). 
In the radio domain, continuum observ ations have been per- 
formed at 20, 6 and 2 cm with the VLA JSandavist et all 19821: 
I Saikia et alJll994rS andavist et al.lll995t) and at 6 and 3 cm by 
Forb es & Norrisl 1119981) with the Australian Compact Array 
(ATCA). These maps show that the nucleus is surrounded by 
several non-thermal radio sources, while the AGN itself is not 
a particularly strong radio source. The radio sources do not co- 
incide with the bright optical knots. A set of NIR images of 
the nuclear region of NGC1365 at very good spatial resolution 
(~0.3") were taken during the co mmissioning of ISAAC at the 
VLT bv lMoorwood et alJ Jl998l) : very red (K-L * 4) counter- 
parts to the brightest radio sources were detected, suggesting 
the presence of deeply embedded sources. 
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NGC1365: radio, MIR and optical 
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Fig. 7. NGC1365: from left to right: radio, MIR and optical images. The map at 3 cm is from iForbes & Norrisl dl998l) : our 
TIMMI 2 map at 11.9 yum is displayed with a non-linear flux scale in order to highlight faint structures and the HST image 
(F814W) is from the archive. On the MIR and optical maps, the locations of the main 3 cm sources are marked by circles of 
radius 0.8". 



Very few MIR obs ervat i ons o f NGC1365 have been per- 
formed to date. lFrogel et al.ldl982l) made a series of MIR mea- 
surements that are only consistent within a factor of two with 
our recent measurements. Such a difference can have various 
origins (differences in filter shape, aperture, and even, over this 
20 years period, genuine variations of the MIR emission of the 
AGN). The map at 11. 9 /mi reveals a 4" extension centred on 
the nucleus (Ml) along the E-W direction (Fig.|2}. This exten- 
sion does not clearly appear on the other maps. At 10.4, 11.9 
and 12.9/mi, the nuclear source (Ml) itself presents a slight 
(2") extension in the E-W direction which is absent in the other 
MIR sources in the field. Therefore, we believe that this exten- 
sion is real. Indeed, to the East of the nucleus, a radio "jet-like" 
structure betw een the nucleu s and t he radio source F has been 
observed bv ISandavist et all i 19951) . Thus the MIR and radio 
emissions of this extended feature might be related. Figure [7] 
compares, at the same scale, the maps in the visible, in the radio 
and in the MIR. Both sources M2 and M3, are clearly detected 
at 10.4 /mi, while at 1 1.9 /mi only M3 shows up and neither of 
them is detected at 12.9 /mi. This indicates that M2 and M3 are 
not strong dust emitters, with a rather flat or decreasing MIR 
spectrum and no PAH emission. The comparison between the 
radio and the MIR maps shows that, as in the case of NGC 1 808, 
the bright radio sources have a counterpart in the MIR, while 
they do not coincide with any strong optical feature. The co- 
incidence between M6 and the radio source G is not perfect, 
if using the map by Forbes & Norris ( 1998). Comparison with 
the 2 cm map bv lSaikia et al.l (119941) gives a much better coin- 
cidence. This kind of discrepancies may come from the uncer- 
tainty in the reconstruction of radio maps when the emission is 
from an extended component together with point-like sources. 
We conclude that sources M4, M5 and M6 are the respective 
counterparts of radio sources D, E and G. The optical HST 
image clearly shows the presence of a dust lane to the north 
of the AGN. At the location of M6, no optical source is de- 
tected, which is understandable as M6 lies in the middle of the 



dust lane. Two optical sources are present at the location of M5 
and might be its counterpart. An extended cone-shaped optical 
source appears at the location of M4, and may be related to it. 
The nature of the MIR/radio sources in NGC 1365 is discussed, 
together with those observed in NGC 1808, in Sect.|3] 



3.6. NGC7469 

NGC7469 exhibits a type 1 AGN hosted in an SBa galaxy. A 
significant fraction of the nuclear luminosity originates in a 
circumnuclear ring of star formation (1 kpc diameter), which 
is seen face-on and thus appears more or less symmetrical 
around the AGN. The circumnucle ar ring has been de tected 
in the optical (speckle imaging bv iMauder et alJll994l) . NIR 
jMarco& Alioinlll998l) and MIR jMiles et alJ 1 1994 1 while 
PAH emission has be en observed by Mazzarel la et alJ (11994) 
and M iles etaljfl994 . Non-thermal radio continuum emission 
is detected both as an extended componen t and as several knots 
located 1.5" to 2" away from the nucle us llUlvestad et alJl 98 1 ; 
ICondon et alJl982tlwilson et alJl99lh . hence belonging to the 
star formation ring. 

The 10.4 /mi flux measured from the TIMMI 2 dataset for 
the central sou rce is in good agre ement with the N-band mea- 
surement from lGoriian et al] ll2004h . The morphology seen on 
the TIMMI 2 map at 11.9 /mi is consistent with high resolution 
NIR images, particularl y the 0.4" reso lution J and H SHARP 
maps published by Gen zel et al.Nl995l) : the circumnuclear ring 
exhibits enhanced brightness to the NE and SW. Surprisingly, 
the 10.4 //m image has a quite different aspect and reveals only 
an extension along the ~N-S direction. A simple interpretation 
is that, except along the N-S axis, the extended MIR emis- 
sion is dominated by PDR emission, with low continuum at 
10.4 /mi and intense PAH emission in the 1 1 .9 //m filter. An ad- 
ditional N-S micro-bar with large dust content could explain 
the observed intense N-S emission at 10.4 /mi. 



10 



Galliano et al.: MIR imaging of active galaxies 



3. 7. Unresolved galaxies 

NGC2992 is an almost edge-on Sa galaxy, interacting with 
NGC2993. It is at a distance of 30Mpc, which gives a scale 
of 150pc/". It harbours a type 2 AGN, located behind a 
prominent dust lane extending along the major axis of the 
galaxy. The radio source at the location of the AGN is pecu- 
liar in shape (shape o f a 'figure 8') and extends up to 450 pc 
JChapman et al.ll2000h . The nuclear region of NGC2992 was 
observed through the TIMMI 2 10.4 and 12.9 fim filters. We did 
not detect any conspicuous diffuse/extended emission down to 
a level of 35 mJy arcsec~ 2 (l<x) at 12.9 fim. The AGN appears 
as a point-like source with FWHM=0.9" on the 10.4 /mi im- 
age. The flu x measurements are in agreement with the mea- 
surements by Goriian et al ] fe004 .The 12.9/10.4 flux ratio has 
a value of 2.44, compatible with both the AGN and PDR tem- 
plates, and thus PAH bands may be pres ent. A detectio n of 
the 11.3 fim PAH band was claimed bv lRoche et ail Jl99ll) . al- 
though Mizuta nT et alJ J 19941) did not confirm this result as they 
did not find the PAH signature at 3.3 fim. 

NGC5995 hosts a type 2 AGN and is located at 100 Mpc, 
leading to a scale of 490 pc/". It has been observed with 
TIMMI 2 through the 8.9 and 1 1.9 fim filters. Comparable flux 
values were found in both filters, suggesting emission from a 
PDR. The 11.9 fim flux is found to be slightly lower than, but 
consistent with, the IRAS flux value, implying that most of the 
MIR emission is from the central source. 

IZwl is the archetypal narrow-line type 1 AGN. It is con- 
sidered as the closest QSO because of its large intrinsic op- 
tical lumi nosity. At a distance of 245 Mpc, the linear scale is 
1200pc/". Schinne rer et alJ ll 19981) find a strong starburst in a 
circumnuclear ring of -1.5" in diameter (1.8 kpc). IZwl was 
observed with TIMMI 2 at 10.4, 11.9 and 12.9 fim. The fluxes 
are found to be in good agreement with the IRAS ones, indi- 
cating that the bulk of the MIR emission comes from the unre- 
solved nucleus. Interestingly, the 12.9/zm flux is significantly 
lower than the 11.9 fim flux (as in the case of PDR emission), 
and is at the same level as the 10.4/im, which is not expected 
for a PDR emission. This suggests the presence of the silicate 
feature in emission. Spectroscopy in the MIR should allow to 
confirm this prediction. 

IIZwl36 is a type 1 AGN, 252 Mpc away (scale of 
1200pc"). It has been observed through the 10.4 and 
11.9 fim filters: equal flux values are found in both filters. Such 
a flux ratio, 1+0.3, is consistent with either the HII region tem- 
plate or the AGN template if the 9.7 fim silicate feature is not 
as deep as in the template. 

4. Active galactic nuclei: building the SED 

This short section focuses on the SED of the AGN in our sam- 
ple. Historically, diagnostics of AGN SED were preferably ap- 
plied to high luminosity objects, in order to avoid complication 
due to the mixing of light from the AGN and from the host 
galaxy. Yet, the study of weaker close-by AGN is important 
since their proximity allows their inner structure to be unveiled 
down to physical scales of a few parsecs, corresponding to the 
expected sizes of some of their components. In such cases, the 



main difficulty arises from the fact that the contrast between the 
AGN and its host galaxy is low, particularly when the AGN is 
surrounded by regions of starburst activity. Then, large aperture 
photometric measurements can greatly differ from the actual 
AGN fluxes. 

As the observational c onstra int s on the toru s mod- 
els jKrolik & Beeelmar] Il986l IPier & Krolikl Il993t 
Granato & Danese 1994; lEfstafhiou & Rowan-Robinsonl 
1995c Kiranato et alJll997l) come from the shape of the SED 
between 1 and 20 /mi, and from the depth of the 9.7 /mi silicate 
feature, it is particularly important to isolate the AGN and 
measure its flux free from contamination by other components. 
We provide a compilation of small aperture measurements 
of the NIR-MIR fluxes for our sample AGN (Table QJ. The 
corresponding SED are shown in Fig. [8] On the top-left panel, 
along with the SED of Circinus, typical SED models have been 
overplotted for type 1 and type 2 AGN, as pred icted using the 
radiative transfer code bv lGranato et alJ l[l997). Type 1 AGN 
show a flat SED (in vL v ), while type 2 AGN show a SED with 
a steep fall off at NIR wavelengths. 

- For Circinus, the compiled SED is from high resolution 
NACO NIR data and our MIR measurements. The SED is 
fully consistent with the ISO spectrum (overplotted) and 
shows that the bulk of the NIR-MIR emission comes from 
a type 2 AGN. 

- For NGC180 8, we plotte d new NIR small aperture mea- 
surements ( Galliano 2005) together with the TIMMI 2 data. 
The MIR points are distributed similarly to those of 
Circinus, although slightly shifted towards lower luminosi- 
ties. This supports the idea that the nucleus of NGC1808 
is indeed an AGN. However, at shorter wavelengths, the 
measurements show a SED increasing towards short wave- 
lengths, suggesting the presence of a bright, stellar core. 
The overplotted la r ge ap erture ISO spectrum is from 
Siebenmorgen et aD J200ll) . Its flux level, about 10 times 
larger than that of the AGN, demonstrates that the central 
1 kpc region is dominated by the emission from the star- 
burst. 

- The AGN in NGC1365 is thought to be a type 2. Hence, 
its SED is expected to be decreasing from the MIR to 
the N IR. New small-aperture NIR measurements ( Galliano 
2005) are plotted together with the TIMMI 2 fluxes and 
show that, in contrary to expectations, the SED is flat and 
consistent with that of a type 1 AGN. The presence of a lo- 
cal maximum at 3.5 /mi suggests a contribution from PAH. 
The AGN in NGC1365 is thus a candidate for the detection 
of PAH around AGN. 

- The NIR-MIR SED of NGC2992 and NGC7469 nuclei are 
respectively consistent with type 2 (steep SED) and type 
1 (flat SED) AGN. 

The remaining sources shown in Fig.[8](NGC5995, IZwl 
and IIZwl36) are described in Sect. 13.71 Detailed and unam- 
biguous modelling of the SE D of these AGN requi res the use 
of additional data: as shown in iGalliano et ai] l2003) photomet- 
ric measurements alone allow the fit of many different models 
of tori. A way to lift this degeneracy would be to use spatially 
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Fig. 8. NIR-MIR spectral energy distributions for the AGN in our sample. Intrinsic luminosities have been calculated using the 
distances reported in Tableland they are provided in Table|3 The vertical bar attached to each point represents the photometric 
uncertainty, while the horizontal bar represents the width of the filter through which the measurement has been performed. Large 
aperture ISO-SWS spectra (thin continuous line) are overplotted for Circinus and NGC180 8. On the top-left panel, we have 
overplotted with thick lines modelled typical SED for type 1 and type 2 AGN (Granato et alJI 19971 Fig. 1, second model). For 
the sake of clarity, the flux level of these model SED has been set at a large value. 



resolved information, which means to observe these AGN with 
IR interferometers. 

5. Circum-AGN environments: embedded star 
clusters in NGC1808 and NGC1365 

MIR imaging at high spatial resolution provides the oppor- 
tunity of detecting young massive dust embedded clusters in 
starbursting galaxies. These clusters emit a bright IR contin- 
uum, infrared nebular emission lines and are strong cm radio 
sources. The IR continuum arises from the dust heated by stars 
in the cluster, while the lines and usually most of the cm contin- 
uum are produced by the HII regions surrounding the massive 
stars in the cluster and/or by supernova remnants (SNR). 

To date, only a very limited number of such embed- 
ded clusters are known, namely the clusters in He2-10 
jKohulnickv & .Tohnsonl fl999t iBeck et al I EoOll IVacca et al J 
2002 |, NGC5253 JTurneret alJ U99 8: Gorii anet all 12001; 
Van zi & Sauyagd l2004 . the Antennae jMirabel et alJ [l998; 
Gilb ert et al.ll2000t). IIZw40 llBeck et alJl2002l) and SBS0335- 
052 JPlante & Saiivagel [2002). The embedded clusters are 
found to be young (a few Myr), massive (~ 10 6 M Q ) and highly 
obscured (A v =5 to 100). 

The fact that the bright compact radio sources surrounding 
the AGN in NGC1808 and NGC1365 are found to be also in- 



tense MIR sources suggests that they are related to very dusty 
objects. Since these sources have no comparably intense optical 
counterpart (except NGC1365/M5), they are very likely to be 
embedded star clusters. Still, this interpretation needs quantita- 
tive support: in the following, we present comparisons between 
the predictions from a simple model for embedded star clusters 
and our observations in the MIR and the radio. We deliberately 
chose to consider a simple model, given the still limited amount 
of data available for comparison. The calculations in this sec- 
tion only intend to give support to the interpretation in terms 
of embedded clusters and to give rough estimates of their prop- 
erties (order of magnitude for their age and mass). Detailed 
and more realistic modelling of such objects is postponed until 
new data, especially spectroscopic measurements, will become 
available. Throughout, the radio spectral indexes are noted a, 
with S Y oc v ". 

5.1. Emission from young star clusters 

Prior to being identified in the MIR, the radio sources surround- 
ing the nuclei in NGC1808 and NGC1365 were commonly in- 
terpreted as radio supernovae or as supernova remnants (SNR). 
This was suggested by the negative values of their cm spec- 
tral indexes, indicating non-thermal emission (the radio mea- 
surements are given in Table |4j. However, the interpretation in 
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Fig. 9. Evolution of the properties of a star cluster. We consider a Starburst99 model for a 10 6 M Q cluster with instantaneous star 
formation and Salpeter IMF. On the plots to the left, evolution of the star cluster intrinsic properties are shown. On the plots to the 
right, predicted observational properties as a function of time are shown. An index a = -0.9 is used for the SNR radio emission, 
(a) Rate of supernova events per year, (b) Non-thermal luminosity at 6 cm produced by the SNR. (c) Lyman photon production 
rate, (d) Thermal luminosity at 6 cm produced in the HII regions, (d) Bolometric luminosity of the cluster, (e) Evolution of the 
10.4 fim luminosity if all the energy of the stars is reprocessed in a black-body at respectively 300 K and 150 K. 



terms of isolated (i.e. not associated with a star cluster) radio 
supernova or SNR can now be discarded since these objects 
would not show up as IR sources for more than a few years 
llWeileretalJl 986) and therefore would have a low probability 
of being detected as such. As we do detect them in the MIR, 
we consider instead that they are radio supernova or SNR as- 
sociated with a young star cluster, where their probability of 
occurrence is large. 

The observed radio spectral index for the sources range 
from steep to rather flat, but always indicate an important 
share of non-thermal processes. Two types of radio emission 
are expected from young compact clusters: free-free thermal 
emission from the HII regions surrounding the massive stars 



(a 0.1 where a can become positive in the case of optically 

thick emission), and non-thermal emission from both radio su- 

pernovae and SNR (a 0.9). We examine hereafter the three 

possibilities. For simplification, we only consider optically thin 
emission, being conscious that such an assumption might not 
always hold true. 



The first mechanism for the production of radio continuum 
is free-free emission in the HII regions surrounding massive 
stars. The observed flux depends on the r ate of emission of H 
ionising photons (E > 13.6eV). Following Turn er et"a"fl JT998 ). 
in the case of optically thin radio emission, the cm flux is re- 
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lated to the number of ionising photons by the following rela- 
tion: 

fel— "(£f(<sns) <" 

where S v is the flux density at frequency v, An is the rate 
of emission of ionising photons, and d is the distance to the 
source. 

Some synchrotron sources appear and disappear within pe- 
riods of months to years following Type II or Type lb super- 
novae events and can be very intense. However, because of their 
short lifetime, and because we observe simultaneously several 
sources in NGC1365 and NGC1808, these radio supernovae 
are unlikely to be the dominant contributors. This contribution 
is therefore discarded in the following discussion. 

Finally, the radio emission from SNR starts about 50 years 
after the supernova event and persists for hundreds to thousands 
of yea rs. We adopt the simple recipe given by Con don & Yinl 
(1990), relating the non-thermal radio emission from galaxies 
to the type II supernova rate: 

(M = lo6 (J_f (jl) (*2l) (2) 

\mJyj SNR \MpcJ \GHz) \yr 'j 

where S v is the flux density at frequency v, a is the spectral 
index, d is the distance of the object and hsn is the supernova 
rate in yr~ ] . This mechanism is bound to play a role and will be 
considered in the discussion. 

Let us evaluate the relative share of the two contribut- 
ing mechanisms, free-free emission and non-thermal emission 
from SNR, along the evolution of a star cluster. As argued 
previously, the observed cm spectral indexes indicate that the 
emission has an important share of non-thermal emission pro- 
duced by SNR. Using the stellar population evolu tionary model 
Starburst99 developed bv lLeitherer etalJJl999l) . we can com- 
pute, for a given cluster, the time evolution of the supernova 
rate and of the emission rate of ionising photons. We choose 
a stellar population model with instantaneous star formation, 
solar metallicity, a total mass of 1O 6 M and a Salpeter initial 
mass function. 

Instantaneous star formation is a straightforward assump- 
tion to model the star cluster evolution. Is solar metallicity 
a reasonable assumption? In spiral galaxies, heavy elements 
are often found to be overabundant -with respect to solar- in 
the central region of the galaxy and also to exhibit an abun- 
dance gradient across the galaxy disc. In the case of NGC1365, 
the oxygen abundance gradient (derived from the analysis of 
HII regions) is rather mild (0.02 dex/kpc) and gives an ex- 
trapolated oxygen a bundance at the centre: 12+log((D/H)= 9.1 
jAlloin et alJll98ll). With a solar oxyg en abundance of 8.66 
in the same unit iAsplund et al.l 120041) . the metallicity at the 
centre of NGC1365 is found to be around 3 times sola r. The 
Starburst99 models presented in iLeitherer et alJ i 19991) show 
that the number of ionising stars and the supernova rate de- 
pend very little on the metallicity at least between Z=0.008 and 
Z=0.04 (solar metallicity being Z=0.02). Hence, our assump- 
tion of solar metallicity for the clusters is well founded. 



The lower and upper mass limits in the IMF are respectively 
set to 1 and 100 M , and supernovae occur for star masses 
>8M . Fig.|5]plots the Starburst99 predictions for the super- 
nova rate (a) and emission rate of ionising photons (c). Fig.[5J> 
and Fig.|5Jl show the non-thermal 6 cm luminosity and the ther- 
mal 6 cm luminosity, computed using formulae and Q, re- 
spectively. The considered mass for the cluster (1O 6 M ) was 
chosen so as to predict radio luminosities which are of the order 
of those observed in the sources of NGC1365 and NGC1808 
(see Table@J. Figs.|5j3 and d show that after 3.5 Myr the non- 
thermal emission overcomes the thermal emission. This epoch 
depends mainly on the evolution of the thermal flux, which in 
turn depends on the slope of the IMF. For the chosen IMF, 
supernovae events stop at about ~ 35 Myr after the birth of 
the cluster. The age of the radio sources in NGC1808 and 
NGC1365 should therefore lie between 3.5 and 35 Myr, given 
their fluxes and spectral indexes. How do the predicted cluster 
properties depend on the IMF cutoffs? 

Decreasing the lower mass limit to 0. 1M would neither af- 
fect the ionising flux of the cluster nor its supernovae rate, but 
the mass of the cluster would be larger by a factor 2.5. On the 
contrary, an increase of the upper mass limit of the IMF will 
have little effect on the cluster masses, but will significantly in- 
crease the ionising luminosity of the cluster at very early stages 
(less than 3 Myr). Yet, even though the supernovae occur ear- 
lier, at these very early stages the cm emission remains dom- 
inated by the thermal emission. Hence, a higher upper mass 
limit for the IMF does not significantly affect the dating of 
the clusters. For the clusters in NCG1808 and NGC1365, the 
cm spectral indexes show that the clusters are at a later evolu- 
tionary stage, more than 3 Myr. Hence, a modification of the 
IMF upper mass would not greatly affect our conclusion. This 
is corroborated by a Starburst99 simulation, made with an up- 
per mass cutoff of 150 M . 

The subsequent step is to check if the observed MIR fluxes 
are consistent with those expected from a star cluster embed- 
ded in dust. Precise modelling would require detailed radiative 
transfer calculations. However, as our aim is to perform a check 
at the level of the order of magnitude only, we can simplify 
the problem and represent the cluster emission as a black-body 
at a given temperature. As mentioned earlier, the continuum 
from VSG observed by ISO in s tarburst galaxi es can be well 
fit by black bodies at T>100K JLaurendll999l) . We consider 
two values bracketing a plausible range of dust temperatures, 
150 K and 300 K. Fig. [9^ plots the bolometric luminosity of 
the model cluster and Fig. |5f translates it into the luminosity 
at 10.4 fim assuming that the whole energy is reprocessed in 
a black-body at 150 K (lower curve) and 300 K (upper curve). 
The 10.4/vm luminosity for T=300K is roughly one order of 
magnitude larger than for T=150K, and in both cases larger 
than the observed 10.4 fim luminosity (Table0Ji. This difference 
is discussed in Sect. 15. 31 

Overall, the interpretation that the MIR/radio sources ob- 
served in the starbursting environments of NGC1808 and 
NGC1365 are embedded young massive star clusters is satis- 
factory. 
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5.2. Dating and weighing the clusters 

Within the framework of the simple cluster modelling dis- 
cussed above, it is possible to derive some quantitative param- 
eters for each cluster. The index of the radio spectrum does not 
depend on the mass of the cluster whereas the supernova rate, 
Lyman photon emission rate and total luminosities are mass- 
scalable quantities. As illustrated in Fig. EH an age can be at- 
tributed to the embedded clusters by comparing the observed a 
value to the predicted one as the cluster evolves. Once an age is 
determined, we can compute the expected total 6cm flux for the 
10 6 M cluster model. Rescaling this flux to the observed one 
provides the cluster mass. After mass and age have been de- 
rived for each cluster, it is possible to predict its 10.4 fim flux 
if, as discussed above, we assume that the total luminosity from 
stars in the cluster is reemitted as a single black-body at tem- 
perature T. Results of these calculations are shown in Table0] 
The final figures we obtain are similar for all the bright 
clusters (except NGC1808/M2). Their ages lie in the 3 to 6 
Myr range. Even though their radio emission is not anymore 
dominated by the HII emission from massive stars, the clus- 
ters are very young. The masses inferred are between 3 10 5 and 
2 1O 6 M , comparable to typical globular cluster masses. 

5.3. Colour-colour diagrams 

Modelling the MIR emission of the clusters with black-bodies 
was an easy and simple way to test this interpretation quan- 
titatively. Yet, as shown with the MIR templates in Sect. 12.41 
emission in the N-band is expected to be more complex, with 
the presence of the silicate feature at 9.7 fim, PAH and emission 
lines ([SIV] 10.4 fim and [Nell] 12.8 /mi). The 10.4 /an filter in- 
cludes the silicate feature and coincides with the [SIV] line, the 
11.9 fim filter coincides with a strong PAH band, and finally 



fa NGC1365 
A NGC1808 



-0.4 



C 




5 6 7 

cluster age [10 s yr] 

Fig. 10. The curve shows the evolution with time of the spectral 
index a^ 6 for our adopted cluster model (see Fig.[9]and corre- 
sponding description in Sect. 15 - 1 1 . For a given model, an age 
may be attributed to every cluster in NGC1808 and NGC1365, 
depending only on the measured 6 value. 



the 12.9//m filter samples the continuum outside the silicate 
feature and coincides with the 12.8 fim [Nell] line. 

Fig.^2p resents two colour-colour diagrams: to the left the 
1 1.9/12.9 flux ratio vs. the 10.4/1 1.9 flux ratio and to the right 
the 10.4/12.9 flux ratio vs. the cm spectral index a$ 6 - 

Let us discuss first the left diagram, comparing the N-band 
colours of the brightest MIR/radio sources in NGC1808 and 
NGC1365 (and the AGN of Circinus) with the N-band colours 
of the template spectra of AGN, HII and PDR fSect.1231. The 
filled symbols feature AGN while the open symbols are asso- 
ciated with the surrounding MIR/radio sources. On the colour- 
colour plot, the MIR/radio sources in NGC1808 cover a narrow 
range of 10.4/11.9 flux ratios. On the contrary, the MIR/radio 
sources in NGC1365 are spread in the 10.4/1 1.9 flux ratio and 
build up a sequence which is offset towards small 11.9/12.9 
flux ratios. This offset could trace the uncertainty in the photo- 
metric calibration of the 12.9 fim image due to local and unreg- 
istered variations of the sky transparency (see Sect. 12. 3> . New 
observations are needed to confirm the 1 1.9/12.9 flux ratios in 
NGC1365. Even considering this uncertainty for NGC1365, 
the MIR colours of the clusters in both galaxies are redder 
than the templates. This could result from the steepness of the 
continuum slope between 10 and 14 /an, and/or from the pres- 
ence of strong emission lines. A steep slope in this wavelength 
range can be produced by a deep 9.7 fim silicate absorption fea- 
ture. The colours of dust screened black-bodies at T=150K 
and 300 K are overplotted in Fig. II lh. Since the 10.4/mi fil- 
ter samples the deepest part of the 9.7 fim silicate feature and 
the 12.9 fim filter samples its outer edge, a large amount of ex- 
tinction (several tens of magnitudes in the V-band) is consis- 
tent with the small 1 1 .9/12.9 and 10.4/1 1 .9 flux ratios observed 
(deep silicate absorption). We are aware that a dust screen 
model is a coarse approximation when the MIR continuum is 
emitted by the dust itself. Yet, it is useful to compare the mea- 
sured 10.4/mi fluxes to those predicted by the model discussed 
in Sect. l5.2l and interpret the differences in terms of extinction 
by a dust screen. The Ay values derived this way range from 
a few magnitudes for T=150K to several tens of magnitudes 
(-50) for T=300K. PAH and/or [Nell] 12.8 fim lines are also 
expected to contribute in producing the large observed ratios. 
Spectroscopic observations will be performed to investigate the 
silicate feature, the continuum slope and contribution of PAH 
and [Nell] lines. 

Fig. II lb provides the 10.4/12.9 flux ratio vs. the cm spectral 
index at 6 . As seen in the previous paragraph, the 10.4/12.9 flux 
ratio decreases as the silicate absorption feature gets deeper, or 
as the [Nell] 12.8 fim line gets more intense. The plot suggests 
a possible correlation between the 10.4/12.9 flux ratio and the 
spectral index. As shown in Sect. 15.21 the spectral index is ex- 
pected to decrease with time. Thus the x-axis of Fig. II lb may 
be regarded as a time indicator increasing from left to right. 
The increase of the 10.4/12.9 flux ratio together with that of 
the cm spectral index is broadly consistent with the picture of 
an evolving cluster, embedded in a cocoon of dust which dissi- 
pates with time. Along with the evolution of the cluster, 

(1) the ionising flux decreases as massive stars disappear 
through supernovae events, leading to a decrease of the [Nell] 
12.8 fim flux, 
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11.9/12.9 vs. 10.4/11.9 10.4/12.9 vs. index a 




0.2 0.4 0.6 0.B -0.2 -0.4 -0.6 -0.8 -1.0 

10.4/11.9 index a 6 

' 3.6 

Fig. 11. (a): The 1 1.9 /urn over 12.9 //m flux ratio vs. the 10.4 yum over 1 1.9 /im flux ratio for the brightest sources in NGC1808 
(triangles), NGC1365 (stars), Circinus (square) and the template spectra (circles). The filled symbol in each object represents 
the AGN. The solid lines show the location of black-bodies with temperatu re T=150K (thin line) and T=300K (thick line) 
with screen extinction ranging from Ay=0 to 100 (extinction law bv lDudlev & Wvnn- Williams! 19971) . The dots along the lines 
represent Ay steps of 10. The last dot to the right corresponds to Ay=100. The error corresponding to our assumed 20% error on 
the photometric measurements leads to errors of 30% for the ratios: the typical error size is shown on the top left part of the plot. 
The small stars show the position that the sources in NGC1365 would occupy if we had overestimated the 12.9 /im by a factor 2 
(see Sect. 12. 31 . (b): The 10.4/12.9 flux ratio vs. the cm spectral index a for the clusters. For each point, the horizontal bar shows 
the +0. 1 error on the spectral index, and the vertical bar shows the 30% error on the flux ratio. The small stars have the same 
meaning as in (a). 



(2) the dusty cocoon dissipates and consequently the dust 
column density decreases. This in turn translates into a weaker 
silicate absorption. 

Both effects result into an increase of the 10.4/12.9 flux ra- 
tio with time. 

6. Conclusion 

New high angular resolution MIR observations have been pre- 
sented for a set of nearby AGN, and one starburst galaxy. Of 
particular interest is a new population of MIR sources discov- 
ered in the central regions of the active galaxies NGC1808 and 
NGC1365 and found to be coincident with cm radio sources. 
These MIR/radio sources are interpreted in terms of young em- 
bedded star clusters. Considering the MIR narrow band colours 
and MIR flux densities, as well as the cm indexes and cm flux 
densities, a first-order model of these sources can be performed. 
From this analysis, we conclude that the MIR/radio sources are 
young (5Myr) globular-like (1O 6 M ) clusters, still embedded 
in their dust cocoon. Further observations and modelling are 
being carried out to better understand their evolution. In par- 
ticular, measurements of the ionised hydrogen line fluxes will 
help assessing the relative contributions of thermal and non- 



thermal radio emission. MIR spectroscopy will provide clues 
about the dust composition. Detailed modelling, including dust 
radiative transfer calculations, is being performed (Galliano & 
Delva, in preparation) and will be used to make new predictions 
about these exciting objects. 
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Table 2. Photometry of resolved MIR sources for each galaxy in our sample in each N-band filter. The symbol '-' indicates that 
the source was not observed through the corresponding filter. Source labelling is according to the scheme described in Sect. [5] 



Target 


Individual 


Flux [mJy]" 


Flux [mJy]" 


Flux[mJy] a 


Flux[mJy] a 


Relative position* 




source name 


8.9 pm 


10.4 pm 


11.9pm 


12.9 pm 


u L J 


u L J 


NGC253 


Ml 






5100 ± 1000 












M2 






400 ± 100 




0.87 


1.04 




M3 






350 ± 100 




1.17 


1.5 




M4 






500 ± 100 




1.99 


1.29 




M5 






250 + 50 




2.73 


2.81 




M6 






250 + 50 




4.15 


3.68 




cr bg 






60 








Circinus 1 ' 


Ml 


7950 + 1600 


6650 ± 1350 


16650 + 3300 


23400 ± 4700 








cr b? 


50 


25 


20 


25 






NGC 1 80S 


Ml 




255 + 50 


520 + 1 20 

\ f — \ 1 _!_ 1 — - V ' 


970 + 190 


n 
u 







M2 




< 10 


40 + 5 


45 + 5 


-1.99 


-3.76 




M3 




10 ±2 


50 ± 10 


80 ± 10 


-4.55 


1.39 




M4 




< 10 


70 ± 10 


170 ± 30 


2.77 


0.54 




M5 




23 + 5 


< 50 


90 + 10 


2.87 


-1.59 




M6 




< 10 


75 ± 10 


145 + 20 


5.82 


-5.25 




M7 




40 ± 5 


125 ± 20 


90 + 20 


5.44 


-7.34 




M8 




15 ± 5 


105 ± 20 


195 ± 30 


2.71 


-5.51 




M9 




30 ± 5 


< 40 


60 ± 10 


0.70 


-3.38 




cr bg 




10 


20 


25 






NGC 1365 


Ml 


410 + 80 


440 + 80 


510 ± 100 


1100 + 220 rf 










M2 


< 40 


20 ±4 


< 20 


< 40 


-4.66 


-5.12 




M3 


< 40 


30 ±6 


30 ±6 


< 25 


-5.40 


-2.56 




M4 


< 40 


< 5 


40 ± 8 


160 ± 30° 


0.41 


7.12 




M5 


< 40 


60 ± 10 


120 ± 20 


320 ± 60° 


2.79 


9.99 




M6 


105 + 20 


50 ± 10 


140 + 20 


505 ± 100° 


4.75 


6.95 




cr bg 


35 


10 


10 


20 






NGC2992 


Ml 




230 + 45 




565 ± 110 








cr bg 




18 




33 






NGC7469 


©=2" 




310 + 60 


565 ± 110 










©=5.6" 


_ 


450 + 90 


980 ± 200 


_ 








crbg 




10 


10 








NGC5995 


Ml 


310 + 60 




300 + 60 










o-bg 


100. 




25. 








IZwl 


Ml 




375 ± 75 


425 + 85 


325 + 65 








o-bg 




10 


18 


16 






IIZwl36 


Ml 




130 + 25 


130 ± 25 










crbg 




13 


13 









None of the individual sources from this table are resolved on the TIMMI 2 images, except tentatively NGC1365/M1 (see Sect. l3.5> . The linear 
upper limits for the source sizes are, given the galaxy distances of TableQ 

NGC253: 12pc; Circinus: 7pc; NGC1808: 40pc; NGC1365: 70pc; NGC2992: 120pc; NGC7469: 250pc for the central peak; NGC5995: 
400 pc; IZwl and IIZwl36: 1 kpc. 

(") Unit for the fluxes is mJy and the unit for cr bg is mJy arcsec~ 2 . For most sources, the quoted error is the 20% calibration error (see Sect. 12.31 
which dominates over the measurement uncertainty on the images. In cases where the measurement uncertainty was not negligible compared 
to the 20% error, it was included. 

(*) The positions are offsets on our MIR images, relative to source Ml, the coordinates of which are set to (0",0"). No significant astrometry 
is possible with this TIMMI 2 data set because of the small field of view. For NGC253, a discussion on the coordinates of Ml (correspond- 
ing to the 0",0" position in this table) is provided in Sect. 13.21 For NGC1 808, Ml can be ide ntified to its counterpart on the cm image 
(see Sect|23J, The B1950 coordinates of the Ml cm counterpart derived bv lSaikia etail <199fj) are RA: 05 h 05 m 58.56 s , Dec:-37°34'36.3", 
with no information concerning the uncertainty. For NGC1365, the B1950 coordinates of Ml radio counterpart are RA: 03 h 31 m 41.80 s , Dec:- 
Se^S^e^' iSandavist et alll995t no information on the uncertainty). This cm image is shown in Sect. 13.51 
( c ) Circinus was also detected at 20 pm, with a flux of 30 ± 6 Jy and cr bg of 1 .2 Jy arcsec" 2 . 

('') The photometric calibration used for NGC 1365 at 12.9 /jm may be greater than the actual one by up to a factor two because of local sky 
transparency variations (see Sect. l2.3> . This would result in flux values that are up to twice smaller. 
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Table 3. Compiled NIR-MIR photometry for the galaxies in our sample. The MIR data presented here are the integrated fluxes 
over the aperture given in the table. 



Name 


Wavelength 


Aperture 


Projected 


Flux 


Reference 


loe(vL ) 




L P 411 J 


[arcsec] 


size [pc] 


rmJvl 




Ters? s — ^1 


Circinus 


1.265 


0.19 


3.6 


1.6±0.2 


Prieto et al. (2004) 


39.9 


2 S Mnc 


1.66 


0.19 


3.6 


4.77+0.5 


// 


40.2 




2.18 


0.19 


3.6 


19±1.9 


// 


40.7 




3.80 


0.19 


3.6 


380±38 


// 


41.8 




4.78 


0.19 


3.6 


1900±190 


// 


42.4 




8.9 


2 


38 


7950±1590 


This work 


42.7 




10.4 


2 


38 


6650+1330 




42.6 




11.9 


2 


38 


16650+3330 


// 


42.9 




12.9 


2 


38 


23400±4680 


// 


43.0 




20 


2 


38 


30000±6000 




42.9 


NGC1808 


2.2 


2 


106 


43.5±4.4 


Galliano (2005) 


41.9 


10.9 Mpc 


3.5 


2 


106 


36.9±7.4 


// 


41.7 




4.8 


2 


106 


19.1+3.8 


// 


41.2 




10.4 


2 


106 


255+51 


This work 


42.0 




11.9 


2 


106 


620+124 




42.3 




12.9 


2 


106 


970+194 


// 


42.5 


NGC1365 


2.2 


2 


106 


78+7.8 


Galliano (2005) 


42.6 


18.6 Mpc 


3.5 


2 


106 


205 ±41 


42.9 




4.8 


2 


106 


177±35 


// 


42.7 




8.9 


2 


180 


410±82 


This work 


42.8 




10.4 


2 


180 


440±88 


// 


42.7 




11.9 


2 


180 


510±102 


// 


42.7 




12.9 


2 


180 


1 10±220 


// 


43.0 


NGC2992 


2.18 


3 


450 


a 2.8±0.2 


Alonso-Herrero et al. (2001) 


41.6 


301. Mpc 


3.8 


3 


450 


8 22.7±2.0 


/; 


42.3 




4.8 


3 


450 


" 35.7±10.7 


// 


42.4 




10.4 


2 


300 


230±46 


This work 


42.9 




12.9 


2 


300 


565+110 




43.2 


NGC7469 


1.1 


3 


960 


8 16.2±2.4 


Alonso-Herrero et al. (2001) 


43.4 


66 Mpc 


1.6 


3 


960 


8 39+5.9 




43.6 




2.22 


3 


960 


a 67.8±17.0 


// 


43.7 




3.80 


3 


960 


" 159±4.8 


Wardet al. (1987) 


43.8 




4.78 


3 


960 


8 259+31 


If 


43.9 




10.4 


2 


640 


310±62 


This work 


43.7 




11.9 


2 


640 


565+113 


// 


43.9 


NGC5995 


2.2 


10 


3200 


38+0.8 


Hunt et al. (1999) 


43.8 




8.9 


2 


980 


310±62 


This work 


44.1 


100 Mpc 


11.9 


2 


980 


300±60 




44.0 


IZwl 


3.5 


8.5 


10000 


1 10+11 


Rieke(1978) 


44.8 


245 Mpc 


10.4 


2 


2370 


375±75 


This work 


44.9 




11.9 


2 


2370 


425±85 


// 


44.9 




12.9 


2 


2370 


325±65 


// 


44.7 


IIZwl36 


10.4 


2 


2440 


130±26 


This work 


44.5 


252 Mpc 


11.9 


2 


2440 


130±26 




44.4 



8 Non-stellar contribution, as computed bv lAlonso-Herrero et alJ fcOOll) 
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Table 4. Measurements and predictions of the properties of the embedded star clusters in NGC1808 and NGC1365. 



Source name 


M2 


NGC1808 
M3 M4 


M8 


M4 


NGC1365 
M5 


M6 


"Observed afj 6 








-0.92 


-0.28 


-0.40 


-0.61 


-0.39 


-0.77 


-0.44 


Estimated Age (Myr) 








>6.3 


3 4 +01 


3 5 +as 


5 - 1+ !'o 


3 - 5+ oi 


5.7 +2 - 5 




"Observed 6cm Flux [10 25 ergs~'Hz 








9.2 


15.2 


14.6 


14.6 


110.5 


58.6 


113.4 


Predicted 6cm Flux [10 25 ergs~'Hz~ 


'] 






9 


57 


64 


41 


64 


37 


53 


b Mass [10 6 Mo] 








1.6 


0.3 


0.3 


0.4 


1.7 


1.6 


2.1 


Observed 10.4/jmLum [10 26 ergs~' 


Hz- 


] 




<15 


15 


<15 


20 


<20 


250 


200 


^Predicted 10.4yumLum [10 26 ergs~ 


l HzT 


'](T= 


=150 K) 


10 


60 


50 


40 


360 


170 


390 


'Predicted 10.4yumLum [10 26 ergs~ 


l Hr 


l lCr= 


=300 K) 


50 


370 


300 


270 


2300 


1100 


2500 



(") The radio measurements from lCollison et all 1 19941) for NGC1808 and lSandavist et alJll995h for NGC1365. The uncertainty on the spectral 
indexes is 0.1. 

(*) Shifting the IMF lower mass cutoff to 0. 1M Q , the mass would 2.5 times larger. 

(') The difference between observed and predicted luminosities at 10.4^m can be accounted for by absorption, the corresponding A v values 
are discussed in Sect. 15,31 



